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' We present a solid-on-solid model of a binary AB compound, where atoms of type A in the top- 

, most layer interact via anisotropic interactions different from those inside the bulk. Depending on 

/"S ■ temperature and particle flux, this model displays surface reconstructions similar to those of (001) 

surfaces of 11- VI semiconductors. We show, that our model qualitatively reproduces many of the 
characteristic features of these materials which have been observed during sublimation and atomic 
layer epitaxy. We predict some previously unknown effects which might be observed experimentally. 
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Apart from a few exceptions (eg. Jl|]), theoretical models of crystal surfaces in nonequilibrium situations have not 
taken into account the effects of surface reconstructions, so far. Consequently, a theoretical understanding of the 
interplay of this phenomenon with the dynamics of growing or sublimating surfaces has not yet been achieved. In this 
paper we will address this important question in the context of a compound material that displays a competition of 

S different vacancy structures in the terminating layer. 
Within the last years, there have been intensive experimental studies concerning the surfaces of II- VI semicon- 
'■^ \ ductors, see eg. |^,[|] for an overview. These studies have revealed the existence of a fairly small number of surface 
reconstructions, which makes them promising candidates for a theoretical modelling. In the following, we will re- 
fer mostly to experimental results for CdTe(OOl) [2-9]. Under vacuum conditions, the CdTe(OOl) surface is metal 
i terminated. The Cd atoms arrange in a vacancy structure, where half of the potential Cd sites are empty. At low 

temperatures, a c(2 x 2) reconstruction is observed, in which Cd atoms arrange in a checkerboard like configuration. 
£NJ , Often a contribution of (2 x 1) arrangement in rows along the (110) direction is observed. At a temperature of 
« 300°C, a reordering occurs, above which a (2 x 1) ordering dominates on a sublimating surface. An additional Cd 
flux stabilizes the c(2 x 2) ordering even at high temperatures. Under a Te flux, the surface is Te-terminated, with 
a (2 x 1) reconstruction. At small Te fluxes, the Te-coverage is 1, at low temperatures and high Te fluxes a higher 
coverage of 1.5 is observed [||j5|J|,p. In this Letter, we present a 2+1-dimensional solid-on-solid model of a binary 
I i compound material to get insight into the following problems: (1) Is the claim Q of an effective thermal equilibrium 
of the surface layer justified under the non-equilibrium conditions of sublimation? (2) What is the effect of an external 
particle flux on the surface reconstruction? (3) How do reconstructions influence techniques of crystal growth like 
■ atomic layer epitaxy (ALE)? 

In the terminating layer of a CdTe(OOl) surface, the potential Cd sites form a regular square lattice, where the 
simultaneous occupation of nearest neighbour (NN) sites in the [010]-direction (denoted as y in the following) is 
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forbidden by electron counting rules |10|. In our model, there is an attractive interaction between nearest neighbours 
in [110](x)-direction (coupling e x ) and diagonal neighbours (coupling Sd). Measuring energy in units of |eJ, we set 
Ed = — 1. A model, where the underlying crystal was assumed to be fixed, has been studied elsewhere flllf . If the 
energy difference between a perfect c(2 x 2) arrangement and a perfect (2 x 1) arrangement is small and positive 
(—2 < e x < —1.9), this model explains the c(2 x 2)-(2 x 1) reordering as an accompanying effect of an order-disorder 
transition in thermal equilibrium. At low temperature T, the system is in an ordered phase, where c(2 x 2) dominates, 
above the transition it is in a globally disordered phase, where the (2 x 1) rows dominate the local environment of a 



H I particle. 

In this paper, we extend this planar lattice gas to a model of a three-dimensional crystal. Currently no values 
of diffusion energy barriers from ah initio or semi-empirical calculations are available, and the existing experimental 
data are insufficient for a parameter fit. Therefore, a realistic modelling of CdTe is beyond the scope of this work, 
which aims at an understanding of fundamental properties of nonequilibrium crystal surfaces. Instead, we consider a 
model of a compound U AB 1 ' 1 with a cubic lattice and a comparatively simple potential energy surface which, although 
not being quantitative, reproduces essential features observed in experiments on CdTc(OOl). Here, the A (B) atoms 
correspond to Cd (Te). Within the framework of a solid-on-solid model, we describe the crystal by a two-dimensional 
array {h xy } x v= i °f integers which denote the height of a column of atoms. We use periodic boundary conditions. 
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The interactions between the particles are described by phenomenological parameters which are supposed to include 
the effects of surface strain and other interactions. To model the layered structure of the CdTe crystal, we assign the 
odd heights to A particles, and the even heights to B particles. Inside the bulk of the crystal, there is an attractive 
interaction between the particles and their NNs which is isotropic in directions parallel to the surface. While there 
is no difference in the interaction of B particles between the bulk and the surface, A particles on the surface interact 
with the anisotropic interactions e x and e d - The other parameters of the model are defined in Figure [l]. Thus, the 
Hamiltonian of the crystal is 

H = e AB n AB + £BB n BB + £bn b AA + e h n h AA + e d n d AA + e x n x AA (1) 

Here, n A s is the number of AB bonds between the layers, ubb is the number of NN pairs of B atoms, n b AA the 
number of A NN pairs inside the bulk, n AA the number of surface A atoms next to a higher column, n AA (n AA ) the 
number of bindings of surface A atoms to diagonal neighbours (NNs in x-direction) at the same height. The stay of 
particles in the wrong sublattice and the formation of NN pairs of surface A atoms in y-direction is forbidden. We 
allow diffusion both to nearest neighbour sites and to diagonal neighbours. Additionally, we permit the diffusion of 
an AB pair with the B atom on top of the A atom, if the diffusion of the B atom alone would end in the wrong 
sublattice. This process is required to preserve the ergodicity of the system. We simulate a Kawasaki dynamics, i.e. 
the rate of diffusion events is vq exp(— (B + AH)/T) if AH > and v exp(— Bq/T) else. The desorption of an A (B) 
particle requires an additional activation energy E A (Eg). This is done using rejection free Monte-Carlo techniques 
employing a binary search tree H] . Particles from an external source arrive at each lattice site with equal probability. 
If the adsorption of the particle at this site would lead to a forbidden state, it is destroyed. 

This dynamics does not depend on the parameter e A b, so it needs not to be specified, and £bb and Eb enter only 
via the sum e' := Ebb +£&• The constant prefactor 1/to := vq exp(— Bq/T) sets the timescale of the model and needs 
not to be specified, if we measure time in units of t . 

In the following, we will use the parameter set e' = e d = —1, e h = 0, e x = —1.97, E A = E B = 1.5 (in units of \e d \)- 
The fact that the difference in the surface energies between perfectly c(2 x 2) and (2 x 1) reconstructed surfaces is 
small is consistent with DFT calculations Q. However, precise values of binding energies as a quantitative input to 
our model are not available. We have verified, that the behaviour of our model remains qualitatively the same for a 
wide range of parameters, as long as the basic property of a small energy difference between the reconstructions is 
preserved. 

To characterize the surface reconstruction quantitatively, we measure the j4-coverage 9 A = n s A rl -/N 2 , and the 
normalized correlations C AA := n\ A /(4n A ur ^') and C AA := n AA /(2n s A r ^'), which quantify the fraction of ^4-atoms on 
the surface which are incorporated in locally c(2 x 2) or (2 x 1) reconstructed areas, respectively. Additionally, we 
calculate LEED intensities to model diffraction experiments. Neglecting multiple scattering, this is done by calculating 
the fourier transform of the array {exp(i Ak hxy)}^ y=1 , where equal diffraction intensities for A and B atoms have 
been assumed. The antiphase condition between different layers is Ak — tt/2. 

Figure ||a shows the evolution of an initially ^-terminated, flat surface under vacuum at a temperature T = 0.5. 
Initially, 9 A increases from zero to an asymptotic value close to 0.5. We find, that this increase can be fitted with 
an exponential relaxation, 9 A = 6 A (1 — exp(— t/r)). Here, r is the time constant of the decay of the 5-terminated 
surface. We find, that the temperature dependence of the corresponding desorption probability p a 1/t follows an 
Arrhcnius law, p oc exp(—E act ,/T), with an activation energy E act . = -Bo+3.0±0.1 After this onset, A and B evaporate 
stoichiometrically, and sublimation proceeds in layer-by-layer mode, which reflects in the oszillations of the specular 
LEED intensity I spec in antiphase (dotted curve in Figure ga). The dependence of the sublimation rate on T follows 
an Arrhenius law with a higher activation energy B$ + 6.0 ± 0.2. Thus, our model reproduces the observation of Cd- 
terminated CdTe surfaces under vacuum. Studies of surface lifetimes |s| and QMS measurements of evaporation rates 
j| have also revealed significantly smaller activation energies for the decay of the Te-terminated surface (~ 0.7eV) 
than for congruent evaporation (1.9AeV in layer-by-layer mode). We note however, that in our model the microscopical 
energy barriers for desorption of A and B are equal. Therefore the difference in the macroscopic activation energies 
is solely due to the stabilizing effect of the surface reconstruction. 

Surprisingly, C AA and C AA oszillate during sublimation, which proceeds layer by layer. Each time a complete layer 
has desorbed and the surface is atomically flat (maxima of I spe c) the c(2 x 2) reconstructed fraction of the surface is 
maximal. On the contrary, a rough surface with a large number of islands (minimal I spec ) seems to prefer the (2x1) 
reconstruction. This can be understood from the fact that the attractive lattice gas interactions are present only 
between particles in the same layer. Thus, the island edges impose open boundary conditions to the lattice gas of A 
atoms on the island. In contrast to a c(2 x 2) reconstructed domain, a (2 x 1) terminated island can reduce the energy 
of its boundary by elongating in x-direction. Since the ground state energies of both structures are nearly degenerate, 
the formation of (2 x 1) may reduce the surface free energy. This picture is confirmed by the fact that islands on 
sublimating surfaces are indeed elongated (See Figure f|a,b) 
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Experiments have revealed a high density of steps on CdTe surfaces due to the intrinsic surface morphology M, 
which leads to a dominant contribution of step flow sublimation. To incorporate this effect in our model, we have 
performed simulations of vicinal surfaces We find, that the oszillations in the correlations dissappear for terrace 
widths smaller than rs 50 lattice constants, where step flow sublimation dominates. Then, the correlations become 
stationary apart from statistical fluctuations. Figure |^a shows the T-dependence of 9 a, C aa , C aa and the sublimation 
rate r su b at a step distance of 32\/2 lattice constants. The steps were oriented at an angle of 45° to the x-axis, which 
corresponds to the preferential (lOO)-orientation of steps observed on CdTe(OOl) M. Snapshots of the surfaces are 
shown in Figure ^. Clearly, at T — 0.55 there is a transition from a c(2 x 2) configuration at low temperature to a 
high-temperature regime where the (2x1) ordering dominates and the material is sublimating. This is the conterpart 
of the phase transition observed in our two-dimensional lattice gas model. In the investigated temperature range 9a 
decreases only slightly from A = 0.498 at T = 0.4 to 9 A = 0.41 at T = 0.8. 

To get insight into the behaviour of a material in an MBE environment, it is important to understand how the 
structure of the surface changes when it is exposed to a particle beam. To this end, we have simulated surfaces 
which were exposed to a flux of pure A or B. Figure [|b,c shows the quantities 9a, C aa and C AA as functions of the 
adsorption rate at a temperature T — 0.4, which is far below the c(2 x 2)-(2 x 1) reordering under vacuum, and at 
T = 0.57, which is slightly above the transition. Applying B fluxes, we can regulate 9a to values smaller than those 
close to 0.5 observed under vacuum. We find, that a decrease of 9a leads to a significant increase in C AA (Fig. ||b,c). 
This effect is analogous to the higher fraction of (2 x 1) ordering on surfaces with a large number of islands: at low 
9 A , ^4-covered islands form on the surface, whose boundaries favour the arrangement of A atoms in rows. On the 
contrary, an A flux increases 9 A and C AA and leads to the formation of a dominant c(2 x 2) reconstruction even at 
temperatures above the reordering under vacuum. The reappearance of a c(2 x 2) reconstruction under a Cd-flux at 
high temperature and the formation of Te-terminated surfaces with coverage 1 even at low Te fluxes are well known 
from experiments on CdTe [9. In these experiments however, the formation of Te dimers on the surface was observed, 
an effect which is not present in our simple model assuming isotropic interactions between B atoms. 

It has been argued [||, that the CdTe surface is close to thermal equilibrium during sublimation, since the bulk 
serves as a particle reservoir which defines chemical potentials for both elements. Comparing our solid-on-solid model 
with the two-dimensional anisotropic lattice gas, we calculate the equilibrium values of C AA , C AA in the planar lattice 
gas model at the value of 9 A we measure for AB at a given temperature, which can be done using transfer matrix 
techniques |ll]] . Although the non-equilibrium conditions of sublimation enhance the dominance of C AA over C AA in 
the high temperature regime, we find qualitative agreement in the behaviour of C AA and C AA . However, a simple 
mapping of the surface layer of the sublimating system on a 2D lattice gas in thermal equilibrium is not possible. 

Atomic layer epitaxy (ALE) provides a scenario for the investigation of our model of the II- VI semiconductor surface 
in a situation present in technical applications. The idea is to obtain self-regulated growth by alternate deposition 
of pure A and B. In the absence of reconstructions, one would expect the formation of complete monolayers of A 
and B during deposition of the elements, yielding growth at a speed of one monolayer (ML) per cycle. However, this 
does not necessarily apply in the presence of surfaces terminated by vacancy structures with submonolayer coverage. 
Experimentally, in CdTe one finds growth rates of w 1 ML/cycle only at T < 260°C, and growth at a speed of « 0.5 
ML/cycle at higher temperatures [g|J|]. Figure ||b shows the evolution of 9 a, C aa and C AA during ALE at T = 0.3. 
Each cycle consists of two phases of length t cyc i = 4 • 10 3 io- In the first phase, an A flux of 5 • 10 -3 ML/£o is applied, in 
the second phase B is deposited at the same rate. Growth was started from a B terminated flat surface. During the 
first A phase the surface becomes ^-terminated with 9 A « 1/2. Since T is far below the transition, the reconstruction 
is c(2 x 2). At the onset of the fi-phase, 9a decreases rapidly. This leads to an increase of C AA , as reported above for 
surfaces under a small stationary 5-flux. At the end of the cycle, «50% of the surface are covered with S-terminated 
islands, since only half of the A atoms needed for a closed monolayer is present. The following yl-phase deposits half 
a monolayer of A again. However, now the reconstruction is preferentially (2 x 1). This is analogous to the maximum 
of C AA we observe in layer-by-layer sublimation after half a monolayer has dcsorbcd, which is due to the influence of 
the island edges on the reconstruction. During the following B phase, A atoms and AB pairs on top of the islands 
diffuse into the gaps between the islands, which leads to the formation of a closed 5-terminated monolayer. In our 
simulations, we observe 2-3 repetitions of this cycle |p^| . There are two effects, which keep it from repeating infinitely. 
First, we obtain a growth rate of 0.44 ML/cycle, which is smaller than the ideal value of 0.5 ML/cycle. Second, 
the diffusion of particles from the islands into the gaps is not complete. Both effects hinder a perfect closure of the 
monolayer at the end of the even cycles, which leads to a damping of the characteristic behaviour. These observations 
closely resemble the ideas developed to explain CdTe growth at a speed of w0.5 ML/cycle in ||. 

The experimental observation of growth rates of 1 ML/cycle in ALE at low temperatures is inconsistent with our 
model in its present form, since the hardcore repulsion between the A atoms hinders the deposition of more than 0.5 
ML of A. However, this effect might be explained with the existence of weakly bound precursor states, where particles 
could stay close to the surface until the deposition of the opposite species allows for an incorporation into the crystal. 
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At higher temperatures this reservoir would be inactive due to the weakness of the binding. Future research might 
also address the influence of an interlayer diffusion barrier on ALE and MBE in our model. Additionally, one further 
step towards a more realistic modelling of CdTe should be the simulation of a zinc-blende lattice which, however, 
requires the specification of a greater number of parameters. 

In summary, we have presented a simple model of a binary compound material which reproduces many effects 
observed in experiments on CdTe(OOl) qualitatively. We conclude suggesting some experiments to test our model. (I) 
The preparation of CdTe(OOl) surfaces with an extremely low intrinsic step density might allow for the observation 
of oscillations in the correlations during layer-by-layer sublimation. (2) The deposition of small amounts of Te on a 
CdTe(OOI) surface should lead to a preferential arrangement of the remaining Cd atoms in a (2 x 1) ordering. (3) 
During Cd deposition in ALE of CdTe a c(2 x 2) reconstruction is expected in cycles starting from flat Te-terminated 
surfaces, while in cycles starting from rough surfaces (2 x I) should be present. 
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FIG. 1. Sketch of a surface showing the interactions between the particles. A particles are shown as dark cubes, B particles 
in light gray. 



4 




FIG. 2. Panel (a): Sublimation starting from a flat B-terminated surface at T — 0.5 and a system size N = 128. Solid lines: 
C d AA , dashed lines: C AA , dotted: I spec . (arbitrary units). For clarity of plotting, the last two curves have been smoothed; C AA 
shows the natural fluctuations. Each oszillation period corresponds to the desorption of one ML. Panel(b): Results of ALE 
simulations: C d AA (solid), C AA (dashed) and the thickness of the deposited film (dotted) at T = 0.3, t cyk i — 4 • 10 3 to and a 
particle flux of 20 ML/cycle. The system size was N = 256. 



5 



0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 
T 




0.1 0.05 0.05 0.1 

Ft 

FIG. 3. The stationary state of sublimation from vicinal surfaces with terraces of 32\/2 lattice constants width, oriented 
diagonally, (a): Temperature dependence of 6 A (squares), C AA (circles) , C AA (triangles) and r su b (diamonds) under vacuum, 
(b), (c) The influence of an external particle flux on these quantities at T — 0.4 (b) and T = 0.6 (c). A fluxes are shown from 
left to right, B fluxes from right to left. Note that an A-flux leads to a re-entry into the c(2 x 2) reconstructed configuration at 
temperatures above the reordering in vacuum. Errorbars are on the order of the symbol sizes. The system size was N = 128. 
Each datapoint has been obtained in a single simulation run. 
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FIG. 4. Surface snapshots. The dots show colums terminated with A atoms, while the B atoms have been omitted for clarity. 
Larger dots denote higher columns. Panels (a), (b): Layer-by-layer sublimation. Panel (a) shows sections of 64 x 64 lattice 
constans of the surface in the simulation run of Figure |a at t = 7000£o (maximal I sp ec), panel (b) at t = 1.9 • 10 4 £o (minimal 
I spec)- Panels (c), (d): Sublimation of a vicinal surface in step flow mode. The pictures show sections of surfaces at the end of 
simulation runs shown in Figure pk. Panel (c): T = 0.5 (below the transition). Panel (d): T = 0.6 (above the transition). 
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